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EM  PHENOMENA  IN  THE  ELF  RANGE 


By 

L.  Brock-Nannestad 


PREFACE  TO  THE  SERIES 


The  series  of  three  reports  published  under  the  genera],  heading  of  "Electro¬ 
magnetic  Phenomena  in  the  Extremely  Low  Frequency  Range"  provides  a  final 
summary  of  the  research  in  this  field  carried  out  by  the  SACLANT  ASW  Research 
Centre  between  the  end  of  1959  and  the  end  of  1964. 

The  reports  are  based  mainly  on  a  number  of  Technical  Notes  distributed  only 
within  the  Centre  itself.  Short  reports  of  current  work  have  appeared  regularly 
in  the  Centre's  Quarterly  Status  Reports  and,  on  occasions,  in  briefing  material 
for  various  Scientific  Committees.  However,  the  present  series  is  the  first 
comprehensive  description  of  the  research  carried  out  in  this  field  by  the  Centre. 

Naturally,  as  in  all  scientific  work,  some  work  that  was  started  turned  out  to  be 
less  fruitful  than  expected.  In  addition,  the  departure  of  scientists  studying 
particular  fields  has  brought  those  aspects  of  the  work  to  an  end.  For  example, 
work  on  a  magnetometer  based  on  nuclear  magnetic  resonance,  which  had  shown 
some  promise,  had  to  be  discontinued  for  this  latter  reason.  In  order  to  maintain 
continuity,  a  number  of  such  inconclusive  or  uncompleted  studies  are  omitted 
from  these  reports. 

The  size  of  the  Electromagnetic  Research  Group  did  not  permit  a  broad  attack  on 
all  the  problems  involved,  and,  after  an  initial  research  period,  the  effort  was 
concentrated  on  a  few  aspects.  For  this  final  summary  these  have  been  divided 
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into  three  groups  —  each  reported  separately. 


The  first  report  deals  mainly  with  EM  fields  of  submerged  dipoles  and  finite 
antennas.  A  large  amount  of  work  in  this  field  has  been  done  by  other 
investigators,  and  some  of  the  results  of  this  work  have  been  used  in  this 
theoretical  study.  The  second  report  deals  with  the  EM  background  noise  in 
the  ELF  range  and  the  instrumentation  developed  for  its  measurement,  which 
constituted  the  main  part  of  the  Centre's  work  on  EM  phenomena  in  the  ELF 
range.  The  third  report,  which  is  classified,  deals  with  the  special  applications 
to  ASW. 

The  major  part  of  the  research  was  carried  out  during  a  period  when  similar 
work  was  being  done  in  many  other  places.  The  number  of  published  papers  on 
the  subject  has  increased  rapidly,  doubling  within  only  the  last  three  years. 

To  keep  track  of  these  developments,  a  bibliography  was  prepared  in  1962. 

Since  then,  the  rapid  increase  has  demanded  a  second  edition  (Ref.  1),  which 
was  published  in  1964.  As  practically  all  phenomena  in  the  ELF  range  are 
related  in  one  way  or  another  to  phenomena  outside  this  range,  the  bibliography 
also  includes  papers  on  geomagnetic  disturbances  below  1  cps  and  VLF  phenomena 
above  3,  000  cps.  It  thus  contains  abstracts  of  about  1,  700  papers  and  is  used  as 
a  general  reference  throughout  this  series.  In  specific  cases,  however,  direct 
reference  is  made  to  the  individual  papers  listed  in  each  report's  bibliography. 

During  the  research  period,  the  following  scientists  (listed  alphabetically)  have 
worked  for  longer  or  shorter  periods  in  the  Group. 

L.  Brock-Nannestad  (Denmark) 

Y.  Faroudja  (France) 

K.  Hansen  (Denmark) 

J.  Marthins  (Norway) 
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P.E.  Mijnarends  (Netherlands) 

H.  Rynja  (Netherlands) 

E.  Sirovich  (Italy) 

T.  Strarup  (Denmark) 

G.  Tacconi  (Italy) 

The  present  series  contains  contributions  from  all  the  above,  but,  for  clarity, 
the  individual  contributions  are  not  presented  separately  and  references  are 
only  made  when  the  work  has  been  reported  in  the  open  literature.  The  series 
of  reports  thus  represents  the  work  of  the  whole  group. 

The  author,  who  has  been  Group  Leader  throughout,  wishes  to  acknowledge 
the  collaboration  and  contributions  of  his  colleagues  in  this  research.  The 
main  part  of  the  work  was  carried  out  during  a  period  in  which  the  Centre  was 
under  the  directorship  of  Dr.  E.T.  Booth  and  Dr.  J.  M.  Ide,  and  the  author 
expresses  his  thanks  to  them  for  their  interest  and  support. 
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THE  ELECTROMAGNETIC  FIELD  FROM  SUBMERGED  DIPOLES 

AND 

FINITE  ANTENNAS 


ABSTRACT 


The  report  is  the  first  of  a  series  of  three  that  discusses  various  EM  phenomena 
in  the  ELF  range. 

The  expressions  for  the  various  components  of  the  EM  field  from  submerged 
dipoles  are  presented.  Although  similar  expressions  can  be  found  in  the 
published  literature,  one  of  the  purposes  of  this  report  is  to  present  them 
together  using  the  same  notations  and  measuring  units. 

Applications  of  the  theories  to  special  situations  such  as  discontinuities, 
transients  and  finite  antennas  are  described  and  the  accuracy  of  field  expressions 
is  determined.  Some  propagation  experiments  are  briefly  mentioned. 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


A 

A' 


amp-m 
amp/ m 


d 

s 


E, 


"E" 


ELF 

ea 

"f" 

f 

%> 

% 

H 


h 

l-  Xo 

* 

L 

'Ao 

M 

m 

mho 


complex  transmission  factor 

correction  factor 

ampere  metre 

amperes  per  metre 

water -depth 

skin -depth 

electric  field-strength;  subscript  indicates  the  particular 
component,  (vpm). 

extremely  low  frequencies,  (1-3000  cps), 
free  space  dielectric  constant 
dielectric  constant  of  conducting  medium 
transfer  functions 
frequency,  (cps). 

intrinsic  propagation  constant,  air 

"  "  ",  conducting  medium 

"  "  "  ,  sea  bed 

magnetic  field-strength;  subscript  indicates  the  particular 
component,  (amp/m). 

depth  of  transmitter 

current  ,  (amp) 

length  of  insulated  wire 

length  of  bare  wire 

free -space  wave-length 

wave-length  in  conducting  medium 

dipole  moment,  (amp-m) 

metre 

unit  of  conductivity 
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/*o 

= 

free  space  permeability 

CO 

= 

angular  frequency 

n 

= 

Herz  potential 

R,  R' 

= 

distance  between  antenna  and  point  of  observation 

r 

= 

characteristic  variable  in  A 

r 

o 

= 

radius  of  insulating  platform 

/> 

= 

horizontal  separation  between  transmitter  and  receiver 

0,  Ot 

= 

conductivity  of  lossy  medium 

t 

= 

time 

vpm 

= 

volts  per  metre 

z' 

= 

co-ordinate  along  antenna 

(x,  y,  z) 

= 

cartesian  coordinates 

(y°  ,  <p  ,  z) 

= 

cylindrical  coordinates 

.*) ) 

= 

elliptical  coordinates 

<e, <*,/>.  Xj> 

= 

parameters  describing  relative  position  of  transmitter  and 
receiver. 
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INTRODUCTION 


A  large  number  of  scientists  have  investigated  the  fields  generated  by  dipoles 
embedded  in  a  lossy  medium.  The  present  author's  "Bibliography  of  EM 
Phenomena"  (Ref.  1)  contains  250  references  on  this  subject.  Most  of  the 
published  papers  are  theoretical,  with  more  or  less  elegant  solutions  of  the 
mathematical  problems  involved.  However,  a  number  of  papers  describe 
experiments  that  confirm  the  basic  propagation  laws.  This  present  report 
leans  heavily  on  the  work  of  a  number  of  scientists,  among  whom  J.  R.  Wait 
has  been  the  most  productive.  In  addition,  during  the  earlier  part  of  the 
research,  the  much-cited  work  of  Banos  and  Wesley  (Ref.  2)  was  used  as  a 
basis  of  calculation. 

Even  if  the  basic  laws  for  the  generation  of  fields  and  for  their  propagation 
seem  to  be  well  established,  there  remain  specific  problems.  In  the  bulk  of  the 
theoretical  papers  the  calculations  have  been  simplified  by  assuming  idealised 
conditions,  such  as  a  semi -infinite  conducting  medium.  The  problem  of 
deviations  from  this  ideal  situation  has  been  approached,  for  example,  by  Wait 
in  his  treatment  of  the  coastline  effects  caused  by  more  or  less  rapid  conductivity 
changes  in  the  horizontal  plane. 

It  is  also  a  general  feature  of  most  of  the  published  papers  that  the  results  are 
presented  in  terms  of  non-dimensional  parameters,  and  that  the  mathematical 
operations  are  subject  to  approximations  that  can  restrict  the  validity  of  the 
expression.  This  can  result  in  confusion  when  comparing  different  papers. 

In  this  present  report  only  the  very  near  fields  of  dipoles  are  considered.  In 
these  cases,  the  range  —  the  distance  between  dipole  and  point  of  observation  — 
is  always  very  much  smaller  than  the  wavelength  in  air,  but  it  can  vary  from 
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zero  to  several  times  the  wavelength  in  the  conducting  medium.  Within  such  a 
range  the  influence  of  the  ionosphere  on  the  fields  is  generally  neglected,  and 
it  is  also  assumed  that  the  static  magnetic  field  of  the  earth  has  no  influence 
on  the  fields.  In  general,  the  conducting  medium  is  assumed  to  be  seawater, 
which  has  a  conductivity  of  4  mho/m. 

At  the  end  of  the  report  a  few  experiments  are  mentioned.  These  were  carried 
out  to  gain  familiarity  with  the  whole  technique,  rather  than  to  prove  or  disprove 
a  theoretical  result.  Model  experiments,  using  small  volumes  of  mercury  and 
high  frequencies  instead  of  large  volumes  of  water  and  ELF,  have  been  suggested 
as  substitutes  for  large-scale  experiments.  This  technique  was  tried  in  an 
attempt  to  determine  experimentally  the  distribution  of  the  currents  around  an 
antenna  several  skin-depths  long  and  at  various  distances  from  the  boundary. 
However,  there  were  considerable  experimental  difficulties  and  the  results 
were  unreliable. 

The  effect  of  reflections  from  the  boundaries  of  the  tank  make  model  experiments 
less  attractive  than  they  appear  on  paper,  unless  a  large  effort  in  lay-out  and 
equipment  is  made.  Similarly,  full  scale  experiments  are  subject  to  the 
influence  of  boundaries  such  as  coastlines.  During  the  earlier  part  of  the 
research  programme  these  effects  introduced  confusing  errors  in  the  meas¬ 
urements,  and  some  effort  was  spent  determining  the  sources  of  error.  A  few 
details  about  this  part  of  the  work  are  given. 

In  this  report  emphasis  will  be  put  on  presenting  the  formulas  and  expressions 
needed  for  the  calculation  of  those  components  of  the  EM  field  that  were  of 
major  interest  to  the  whole  project.  As  shown  in  the  second  report  (T.  R.  37) 
of  this  series,  the  horizontal  electric -dipole  is  theoretically  the  most  efficient 
receiver  for  EM  energy  in  the  ELF  range  when  working  from  a  floating  platform, 
and  for  this  reason  theoretical  work  was  also  concentrated  on  the  horizontal 

|^§£tric -dipole  and  the  field  generated  by  it. 

k 
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1. 


THE  EM  FIELD  FROM  DIPOLES 


Many  authors  have  dealt  with  the  problem  of  calculating  the  EM  fields  from 
dipoles.  Probably  Sommerfeld  (Ref.  3)  was  the  first  to  give  the  formal 
solution.  Based  on  this  work.  Wait  (Ref.  4)  gave  self-contained  derivations 
for  the  special  cases  that  arise.  The  following  rfesumfe  will  be  based  mainly 
on  this  latter  work. 

The  field  from  dipoles  is  generally  divided  into  three  overlapping  distance 
ranges  that  relate  the  horizontal  separation,  ^  ,  between  transmitter  and 
receiver  to  the  free-space  wavelength,  A0  : 

Case  I,  />»  ,  is  the  "radiation  field". 

Case  II,  p  x,  ,  is  the  "near  field",  and 

Case  III,  /£>  9\0  ,  is  called  the  "quasistatic  case". 

In  the  present  research  on  EM  phenomena  the  last  case  is  the  most  important 
because  the  free-space  wavelengths  in  the  ELF  range  are  extremely  large.  The 
solutions  in  Case  III  are  also  dependent  on  the  relation  between  the  horizontal 
separation,  ,  and  the  wavelength  in  the  conducting  medium,  \  This 
relation  has  been  treated  in  detail  by  Wait  (Ref.  4).  The  derivations  will  not 
be  repeated  here;  only  the  final  results  and  the  ranges  within  which  they  are 
valid  will  be  presented. 

1.  1  Sinusoidal  Excitation 

In  field  problems  of  this  kind  it  is  customary  to  let  the  conducting  medium  occupy 
the  half -space  z  o  with  reference  to  a  Cartesian  coordinate  system  (x,  y,  z). 
The  electric  dipole  of  infinitesimal  length  d£  is  oriented  along  the  x-axis  and 
situated  at  (o,  o,  h). 
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The  intrinsic  propagation  constant  of  the  upper  medium,  air,  is  defined  as: 


=  -  e0/Mo  co- 

\  i  .  ^  - 


(Eq.  1) 


that  of  the  conducting  medium  as: 


f 


2 

4  ^  -i  0}^  CO  ~  £1/Uc 


CO 


(Eq.  2) 


where  C.Q  and  are  the  free  space  dielectric  constant  and  the  free  space 

permeability,  respectively,  and  where  and  are  the  conducting 

medium  dielectric  constant  and  the  conducting  medium  conductivity,  respectively. 


In  the  frequency  range  of  interest,  the  first  term  in  Eq.  2  is  much  larger  than 
the  second  term;  therefore,  throughout  this  section,  the  second  term  will  be 
neglected.  In  other  words,  this  means  that  the  conduction  current  is  considered 
much  larger  than  the  displacement  current  in  the  conducting  medium.  Eq.  2 
then  reduces  to: 


i  +_i_ 

<5 


JZ 

O' 


+  4- 

e 


(Eq.  3) 


where  Q  is 

c o =  2  rff 

frequencies 


the  skin-depth  in  the  conducting  medium  at  the  angular  frequency 
.  Fig.  1  represents  the  skin-depth  as  a  function  of  the  relevant 
for  various  conductivity  values  typical  of  seawater.  It  is  further 
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assumed  that 


%  »  X 


which  will  be  true  in  the  case  of  seawater. 


For  the  range  of  horizontal  distance,  p  ,  in  which  simultaneously  «  1 

and  i*v>i  »  1  ,  the  expressions  for  the  field -strength  of  the  various 


components  of  the  horizontal  electric -dipole  are  as  given  in  Fig.  2.  This 
situation  is  generally  considered  to  be  the  most  important,  but  it  is  interesting 
to  compare  these  expressions  with  the  corresponding  expressions  for  a 
horizontal  magnetic -dipole, given  on  the  same  figure.  Similarly,  f^r  the  sake 
of  completeness.  Fig.  3  gives  the  expressions  for  both  vertical  electric  and 
vertical  magnetic  dipoles. 

In  all  expressions,  the  subscript  1  has  been  omitted  from  both  the  propagation 
constant  and  the  conductivity  symbols.  Thus  *Y  =  ^  and  C  -  •  The 

source  strength  of  the  dipole  is  expressed  as  the  dipole  moment  M,  which  is 
defined  as: 


M  =  I  •  dl 


for  the  electric  dipole,  and  as 


M  =  I  •  dA 


for  the  magnetic  dipole,  where  dA  represents  an  infinitesimal  area.  Time 
enters  according  to  the  factor  exp(i  oo  t),  but  this  has  been  left  out  in  the 
present  part  of  this  report. 
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The  expressions  for  the  field  of  the  horizontal  electric -dipole  become  more 
complicated  when  the  condition  jvi  =  r vi  »1  is  not  fulfilled.  In  this 
case  only  the  components  E  p  ,  E  tp  and  can  be  written  in  concise 
form.  According  to  Wait  (Ref.  4): 


E. 


v5 


2  7Ttf 


[l  +  (1  *  v)  e  T/!|  COS  <p 


(Eq.  4) 


V  -  JX-  .X 


LT  M  .  _1 

.n2  />4 


2  ~(l  p)e  sin  f 

3  -  (3  +  >  + 'i>2)  e"'^ 


Stn 


(Eq.  5) 


(Eq.  6) 


These  equations  are  valid  for  h  =  z  =  0  and  Ivl  «  1  ,  The  value  of  |^/^| 
is,  however,  unrestricted.  For  1%*!  » i  the  expressions  agree  with  those 
given  in  Fig.  2,  and  when  Ivl  approaches  zero  they  reduce  to  forms 
consistent  with  potential  theory.  The  fields  do  not,  in  this  case,  vary  in  a 
simple  manner  with  z  and  h.  Nor  does  the  exponential  attenuation  with  depth  — 
which  was  present  in  the  expressions  for  M  ^ >  i  —  apply  to  Eqs.  4,  5 
and  6.  Wait  (Ref.  4)  presents  expressions  for  the  fields  in  terms  of  modified 
Bessel  functions  for  the  case  where  z  >  o  and  h  y  o. 


An  attempt  has  been  made  to  verify  Eqs.  4  and  5  experimentally  and  to  determine 
the  depth  dependence.  In  a  full  scale  experiment  one  dipole  is  used  as  the 
transmitter,  and  an  electrode  system  —  of  the  type  described  in  the  second  report 
(T.  R.  37)  of  this  series  —  is  used  to  measure  the  horizontal  field  close  to  the 
surface.  During  work  at  sea  it  can  be  difficult  to  maintain  the  proper  orientation 
for  the  measurement  of  or  To  su^  these  particular  conditions,  Eqs.  4 
and  5  have  been  modified  in  the  following  way. 


12 


Fig.  4  shows  the  geometry  as  seen  from  above  the  surface.  The  transmitting 
dipole  is  placed  at  the  origin  of  the  x,  y  co-ordinates.  The  field-strength  is 
measured  along  an  arbitrary  axis  x^,  that  makes  the  angle  6  with  the  x-axis. 
For  the  purpose  of  calculation  it  is  assumed  that  the  transmitter  is  fixed  and  that 
the  receiver  moves  along  the  x^  axis.  (During  experiments  it  was  found  better 
to  keep  the  receiver  fixed  and  to  move  the  transmitter  along  the  x^  axis,  but  this 
does  not  change  the  calculations. ) 


The  various  distances  and  angles  needed  for  the  calculation  are  shown  in  Fig.  4  , 
and  the  relationships  between  some  of  them  are  given  below.  The  zero  point 
of  the  x^  axis  is  placed  where  the  transmitter  and  receiver  are  closest,  the 
distance  between  them  being  po  .  With  the  symbols  used  in  Fig.  4  one  has: 


Using  Eqs.  4  and  5,  the  field -strength  along  the  x^-axis  is  expressed  as  follows: 


v  -  m  i 

^*1  27fcr  /=>* 


cos  0+«)  “sin^sinOf  4-^  cos  (*»-«)]  (Eq.  7) 


where  A  =  (l+yp)e 


-v 


Introducing  the  above  relations  between  the  angles  into  Eq.  7,  gives: 


rp  _  i  .  _J_ 
r  “  4lftf  P* 


2Acos  0  -  cos  0  —  3 cos  (e-2„ui) 


(Eo.  3) 
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where  F  =  E^/M  is  called  the  transfer -function  and  expresses  the  field -strength 
per  unit  of  dipole  moment.  It  is  expressed  in  the  units: 


vpm 

amp-m 


or 


ohm/m 


2 


According  to  Eq.  8,  F  is  completely  determined  by  the  geometrical  configuration, 
as  given  by  6.  P.  and  Xj.  The  complex  factor  A  in  Eq.  7  depends  both  on  the 
geometry  (  P  )  and  on  the  frequency,  and  has  to  be  evaluated  when  a  numerical 
value  of  F  is  required.  Eq.  8  shows  the  interesting  fact  that  F  becomes  frequency 
independent  for  Q  =  . 

Ct 


Factor  A  can  be  evaluated  as  follows: 


with 


A 

'X 


(1  +  e~T/C> 

i 

8 


and  some  rearrangements 


j 


(Eq.  9) 


The  skin-depth  5  is  expressed  by  Eq.  3, 
in  Eq.  9  can  be  expressed  as- 


so  that  the  characteristic  variable 


p/b 


r 


0.0019869  P 


which  transforms  Eq.  9  to: 


(Eq.  10) 


(Eq.  11) 
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Eq.  11  is  introduced  into  Eq.  8,  from  which  the  transfer  -function  is  calculated. 

A  programme  was  written  for  the  Elliott  503  computer  to  evaluate  the  modulus 
of  the  transfer-function  in  the  range  of -200  —  200  m  for  various 

combinations  of  the  frequency  f,  the  distances  of  closest  approach  Po  ,  the 
conductivity  ($  and  the  angle  Q  .  The  programme  is  reproduced  as  Appendix  1, 
and  Fig.  5  shows  some  examples  of  calculated  transfer -functions. 


Even  if  the  case  of  0  = 


r c 


is  covered  by  Eq.  8,  it  is  worthwhile  to  consider 


this  case  separately.  If  the  field -strength  in  the  y-direction  is  called  E^,  then: 


E 


3M 

27Ccr 


(Eq.  12) 


or,  if  expressed  in  terms  of  and  x^. 


(Eq.  13) 


This  expression  is  comnletely  symmetrical  in  /O0  and  x^,  and  shows  that  E^  is 
zero  when  either  or  Xj  equals  zero. 

It  can  furthermore  be  shown  that  E  ,  has  a  maximum  at  the  point 

yi 

X1  =  "2  P.  (Eq-  I4> 

These  interesting  characteristics  of  this  component— including  its  independence 
of  frequency— would  facilitate  the  interpretation  of  experimentally-determined 
field-strength.  Unfortunately,  the  practical  difficulties  of  keeping  the  orientation 
of  transmitter  and  receiver  at  exactly  90°  are  rather  large. 
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From  an  experimental  point  of  view  it  is  worth  considering  the  relative  magnitude 
of  the  various  horizontal  components  and  their  dependence  on  JQ  .  The 
expressions  from  Fig.  2  have  been  used  for  an  illustration,  as  shown  in  Fig.  6. 
The  relative  magnitude  is  expressed  in  terms  of. 


il 


(A  [ 
2 7\o  '  p* 


(Eq.  15) 


where  "E"  represents  the  components  as  illustrated  in  the  lower  left-hand  corner 
of  Fig.  6.  The  expressions  are  related  in  the  following  way: 

M  _  «  «  _ ,  "  \ 

ip  -  Cos  rp  ,  r<p  =  2  sin  cp  \ 

"Fp  -  'Fp*  =  0-9  jp  <p  =  26.  t>° 


In  connection  with  Eq.  16  and  Fig.  6,  it  is  worth  mentioning  that,  with  two 
perpendicular  electrode  systems,  it  is  not  possible  to  determine  the  direction  from 
the  point  of  observation  to  the  transmitting  dipole  by  means  of  a  single  meas¬ 
urement,  unless  the  relative  orientation  between  the  transmitter  and  the  electrode 
systems  is  already  known.  This  property  of  the  near  fields  creates  the  practical 
problem  of  keeping  track  of  position  and  orientation  during  experiments. 
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1.  2  Corrections  for  Discontinuities 


So  far,  all  the  expressions  given  for  the  field -strength  have  been  computed  for 
the  condition  of  a  conducting  half-space.  Since  this  condition  cannot  always  be 
fulfilled  during  experiments,  two  special  cases  of  discontinuities  have  been 
considered  and  the  appropriate  corrections  found. 

The  first  case  is  the  influence  of  the  sea-bed  on  the  propagation,  and  thus  on  the 
field-strength  as  measured  at  the  point  of  observation.  The  model  considered 
consists  of  a  conducting  sea  of  conductivity  Oj  =  O’  and  of  depth  d.  Fig.  7 
illustrates  the  situation.  The  formal  approach  is  to  derive  the  fields  from  a 
Hertz -vector  having  x  and  z  component  and  satisfying  the  proper  boundary 
conditions  at  z  =  o  and  z  =  d. 

The  rather  complicated  expressions  can  be  simplified  if  one  assumes  that  the 
conductivity  of  the  sea-bed  is  sufficiently  low,  and  that  the  ranges  are  small 
compared  with  the  wavelength  in  air  and  in  the  sea-bed.  With  reference  to  Eqs.  1 
and  2  this  means  that  ^  0  .A  further  simplification  is  made  by 

assuming  that  the  range  is  large  compared  with  the  water  wavelength,  and  that 
the  transmitter  and  the  receiver  are  at  the  surface. 

Under  these  circumstances  it  is  found  that,  to  the  first  order  approximation,  the 
sea-bed  at  depth  d  contributes  a  correction  factor 


(Eq.  17) 


to  the  field  for  an  infinitely  deep  sea.  The  transmission  path  for  this  situation  is 
also  shown  in  Fig.  7.  The  direct  path  from  the  transmitter  to  the  receiver  is  the 
same  as  before.  The  indirect  path,  where  the  signal  travels  through  the  sea-bed, 
contains  a  total  water-path  of  length  2d.  Apart  from  the  phase -shift,  this  part  of 
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the  signal  will  receive  an,  attenuation  of  e 


-2d 


,  where  £  is  the  skin-depth. 


8 

The  factor  2  in  front  of  the  exponential  can  be  explained  by  referring  to  Fig.  7. 

If  the  transmitting  dipole  is  raised  towards  the  surface,  its  first  image,  as 
shown  by  the  broken  line,  is  lowered  towar  ds  the  surface  and  finally  coincides 
with  the  dipole.  Seen  from  the  conducting  medium  this  means  that  the  strength 
of  the  dipole  is  doubled.  (This  point  of  view  agrees  with  optics:  when  a  light 
ray  is  reflected  by  a  denser  medium  there  is  a  phase-shift  of  180°.  However,  if 
the  ray  comes  from  the  denser  medium  and  is  reflected  by  the  interface  with  a 
less  dense  medium  no  phase-shift  occurs. ) 

The  correction  factor  of  Eq.  17  is  written  in  the  following  way: 

A  -  (Eq.  18) 

A  plot  of  the  function  A'  is  shown  in  Fig.  8.  The  modulus  is  given  by  the  length 
of  the  vector  from  the  origin  to  a  point  on  the  curve,  which  is  labelled  with  various 
values  of  a  =  2d/  8  .  As  seen  from  the  curve,  there  is  a  considerable 

enhancement  of  the  signal  when  the  water-depth  is  small  in  comparison  with  the 
skin-depth.  The  factor  has  a  minimum  for  a 2.  5  (d  ^  1 . 25  5  ). 

The  second  case  of  disc  ntinuity,  which  will  be  discussed  briefly,  is  the  influence 
of  the  platform  used  for  the  measurement  of  EM  fields.  During  the  earlier  parts 
of  the  programme,  the  island  of  Tino  in  the  Gulf  of  La  Spezia  was  the  site  of  two 
measuring  stations.  A  more  detailed  description  of  the  equipment  used  for  the 
measurement  of  EM  background  noise  is  given  in  the  second  report  of  this  series 
(T.  R  .  37).  The  stations  were  also  used  during  propagation  experiments  with 
the  same  equipment,  when  one  of  the  horizontal  magnetic  components  in  air  was 
measured.  The  influence  of  the  island  on  the  fields  to  be  measured  has  been 
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calculated  in  terms  of  a  very  simple  model.  The  following  is  a  brief  outline  of 
the  procedure,  together  with  some  final  results. 


The  treatment  of  the  problem  makes  use  of  the  model  shown  in  Fig.  9,  and  follows 

the  procedure  of  Kaden  (Ref.  5).  The  model  consists  of  an  infinitely  thin  plate 

having  infinite  conductivity  and  containing  a  hole  of  radius  r  .  In  the  half -space 

above  the  plate  there  is  a  vertical  electric  field  E  ,  which,  at  great  distance 

z 

from  the  hole,  has  the  value  E  .  The  field  close  to  the  hole  will  be  disturbed 

o 

by  the  presence  of  the  hole,  with  the  result  that  the  field  penetrates  into  the  half¬ 
space  below  the  plate.  The  problem  consists  of  computing  the  field  E  in  the 

z 

neighbourhood  of  the  hole. 

In  relation  to  the  actual  situation  at  the  island  of  Tino  this  crude  model  is 

acceptable.  The  sea  around  Tino  is  10  -  30  m  deep,  which  is  small  compared  with 

the  horizontal  dimensions  of  Tino  —  which  can  be  approximated  by  a  circle  of 

radius  r^  =  200  m.  In  comparison,  the  sea  can  thus  be  considered  as  an  infinitely 

thin  plate.  As  to  conductivity,  it  is  assumed  that  Tino  has  nearly  the  same 

properties  as  air,  so  that  the  conductivity  of  the  sea  is  very  much  greater  than 

that  of  the  island.  Apart  from  the  top  layer,  the  island  consists  of  marble,  the 

-8 

conductivity  of  which  is  about  10  times  that  of  seawater. 


Kaden  (Ref.  5,  p.  222)  gives  the  solution  to  the  problem  as  a  closed  expression 
in  elliptical  co-ordinates: 


E  =  grad  X 


where 


v  -  .  n 

A  '  7T  T[ 


(Eq.  19) 
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in  which  the  elliptical  co-ordinates  (  ?.n  )  are  related  to  the  cylindrical  co¬ 
ordinates  (  &  t  z)  in  the  following  way: 

^  -  r:(l^f)(l-f) 

o  ^  Tj  —  1  ?  —  OO  <  <  OO 

The  (p  -co-ordinate  is  neglected  because  of  the  symmetry  of  rotation. 

The  z-component  of  the  gradient  is  then  expressed  in  terms  of  the  elliptic 
components,  yielding  the  following  final  expression: 


(Eq.  20) 


(Eq.  21) 


In  a  similar  way,  the  horizontal  magnetic  field  parallel  to  the  plate  is  given  by 
Kaden  (Ref.  5,  p.  250)  as 


H.  =  grad  X 

x  =  tH,  p[i  * 


(Eq.  22) 


where  <p  is  the  angle  between  the  radius  vector  and  the  direction  of  the  magnetic 
field  Expressing  the  x-component  in  terms  of  elliptical  components  yields  the 
final  expi  espion 


aocttj." 


1+ 


1_ 

7T 


(Eq.  23) 


which  is  very  similar  to  Eq.  21. 


Numerical  values  for  certain  points  of  observation  have  been  computed  on  the 
basis  of  an  approximation  of  the  island's  proportions.  Fig.  10  shows  the  west- 
east  cross-section  of  Tino  in  full  line  and  indicates  the  positions  of  the  two 
stations,  A  and  B.  The  approximation  used  is  shown  by  the  dimensions  of  the 
dotted  lines. 

Fig.  11  shows  the  corresponding  values  for  H^/H^  and  E^/E^  at  various  points 

on  the  model  of  the  island,  the  co-ordinates  (  ^  ^  )  of  which  are  given  in 

brackets.  The  most  interesting  part  is  H^/Hq,  which  has  a  value  of  about  0.  76 

at  station  A  and  of  about  0.  5  at  station  B.  This  has  been  observed  experimentally 

by  simultaneous  measurements  of  the  magnetic  component  of  the  background 

noise  at  A  and  B.  The  value  of  E  /E  at  station  A  is  close  to  1.  This  approach 

z  o 

does  not  permit  the  calculation  of  the  ratio  at  station  B  and  another  model  would 
have  to  be  considered  for  this  case. 

Similar  effects  occur  where  the  measuring  place  is  surrounded  by  mountains  — 
in  fjords  and  lochs  for  example. 

1.  3  Transient  Excitation 

So  far,  the  fields  have  been  calculated  on  the  basis  of  sinusoidal  excitation  of  the 
dipole.  Many  of  the  phenomena  to  be  studied  —  the  background  noise,  for 
example  —  are,  however,  of  impulsive  character  and  it  was  thought  necessary  to 
study  the  medium  under  the  condition  of  transient  excitation  of  the  dipole. 

The  treatment  was  based  on  Laplace  transform  theory,  and  assumed  a  dipole 
excited  by  a  step -function  current  in  the  conducting  medium.  A  paper  by 
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Mijnarends  (Ref.  6)  gives  details  of  the  treatment  and  the  limits  of  validity  of 

the  expressions.  Besides  the  electric  and  magnetic  components  in  the  conducting 

medium,  the  vertical  electric  field  in  air  E  at  various  points  on  the  surface 

oz 

was  examined.  Since  the  expressions  are  rather  involved,  numerical  results  are 

given  only  for  E  and  H„  . 

&  J  oz  lx 

The  response  of  E  is  a  pulse  v/ith  a  rise -time  proportional  to  the  square  of  the 

oz 

depth  of  the  dipole  but  independent  of  the  horizontal  distance,  and  with  a  long 
trailing  edge.  The  magnetic  component  Ht  rises  slowly  to  a  stationary  value. 
Its  rise-time  increases  with  increasing  horizontal  distance.  It  is  also  shown  that 
there  is  a  marked  difference  in  attenuation  with  depth  —  and  sometimes  also 
with  horizontal  distance  —  between  transient  and  sinusoidal  signals.  Some  of 
the  results  obtained  could  be  checked  against  the  results  obtained  by  Wait  (Ref.  7) 
who  used  exponential  and  bell-shaped  source  functions. 


During  the  course  of  this  study  it  was  also  found  that  the  frequency  response  of 
the  detecting  equipment  has  to  be  taken  into  account  when  determining  the 
range  dependence.  This  is  of  importance  for  propagation  experiments.  The 
following  general  hypothesis,  the  validity  of  which  has  not  been  proven  experimen¬ 
tally  due  to  lack  of  time,  has  been  suggested.  Consider  a  transient  disturbance  of 
the  electromagnetic  field  propagating  over  a  highly  conducting  half -space.  If  it 
depends  on  no  other  distance  than  the  horizontal  distance,  P  ,  between  the 
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transmitter  and  detector,  and  if  its  amplitude  is  proportional  to 


—oC 


,  then, 

/3 


if  it  is  detected  with  a  detector  having  a  frequency  response  proportional  to  OJ 

-  (<x:  t-  2/3) 

the  output  of  the  detector  will  depend  on  horizontal  range  as 


It  was  shown  in  the  previous  section  that  the  field  is  influenced  by  the  presence 
of  the  sea-bed  and  the  non-conducting  island.  It  was  also  shown  that,  in  the 
sinusoidal  case,  the  influence  of  the  sea-bed  could  be  expressed  in  terms  of  a 


<>  r\ 


correction  factor  depending  on  the  ratio  of  the  water-depth  to  the  skin-depth. 

A  similar  influence  will  be  present  in  the  case  of  transient  excitation.  If  the 
transient  is  expressed  in  the  frequency  domain  it  is  seen  that  the  presence  of 
the  sea-bed  will  enhance  the  amplitude  of  the  low-frequency  part  of  the  signal, 
where  the  water-depth  is  small  compared  with  the  skin-depth.  The  high 
frequency  parts  will  be  attenuated  in  the  water -path,  and  there  will  be  no 
change  from  the  half-space  condition.  The  waveshape  at  the  point  of  observation 
will  thus  look  different  from  the  waveshape  computed  on  the  basis  of  Ref.  6. 
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2 .  THE  EM  FIELD  FROM  FINITE  ANTENNAS 

The  radiation  from  a  finite  antenna  in  a  conducting  medium  has  received 

attention  from  several  authors.  In  general,  the  antenna  is  described  as  an 

insulated  wire  connected  to  the  conducting  medium  through  electrodes  that  have 

small  dimensions  but  that  effectively  short-circuit  the  wire.  Under  these 

circumstances  the  dipole  moment  of  the  antenna  is  I  .  £  ,  where  the  current,  I  , 

o  o 

is  constant  along  the  wire  and  Ji  is  the  length. 

The  fields  from  such  an  antenna  can  thus  be  computed  in  the  same  way  as  those 

_  / 

from  the  infinitesimal  dipole,  provided  that  the  horizontal  distance  ,  as 
measured  from  the  nearest  part  of  the  antenna,  fulfills  the  conditions  of  the 
propagation  law  used.  The  contribution  of  the  electrode  to  the  dipole  moment  is 
generally  neglected  by  assuming  that  the  electrodes  are  spheres  having  a  radius 
that  is  small  in  comparison  with  the  skin-depth. 

However,  spherical  electrodes  are  not  used  for  practical  transmitting  antennas 
at  sea.  Instead  one  uses  the  so-called  "short-circuited  coaxial  antenna",  in 
which  the  extension  of  the  centre -conductor  is  in  contact  with  the  sea.  Such  an 
antenna  was  described,  for  example,  by  Moore  (Ref.  8),  who,  however,  also 
neglected  the  radiation  from  the  bare  wire.  This  is  not  generally  per  missible, 
and  an  attempt  was  made  to  determine  this  radiation  and  to  find  the  corresponding 
dipole  moment. 

This  problem  was  approached  according  to  Fig.  12,  which  shows  a  coaxial 
antenna  oriented  along  the  z'-axis,  with  its  insulated  part  extending  from  the 
origin  to  and  the  bare  parts  from  £  to  Al.  The  point  of  observation  is 
P,  and  the  distance  between  P  and  a  point  z'  is  approximated  by: 


t 


(Eq.  24) 
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which  is  valid  for  R  and  R1  )£>  £+  L.  (The  drawing  in  Fig.  12  does  not  fulfill 

this  condition). 

For  an  elementary  dipole  one  has  the  Herz  potential 


I0  l  -VI 

47n?cr  € 


(Eq.  25) 


For  the  insulated  part  of  the  short-circuited  coaxial  antenna  this  gives: 


ins 


Ig  ^ _  /  T€coS0 

4-TTcr'Kycos0  ' 


i) 


(Eq.  26) 


For  the  uninsulated  part,  the  attenuation  of  the  current  is  the  same  as  that  of  a 
plane  wave  in  sea-water.  Since  L  is  rather  small,  however,  there  will  be  a 
reflection  from  the  end  at  z1  =  ■£  +  L;  this  depends  on  the  load  impedance,  which 
is  rather  uncertain  and  is  difficult  to  calculate. 


For  estimating  the  magnitude,  two  extreme  cases  are  considered:  (1)  where  the 
bare  wire  is  treated  as  open-circuited  transmission  line,  and  (2)  where  the 
bare  wire  is  equivalent  to  a  short-circuited  transmission  line.  In  both  cases 
the  current  I(z')  in  the  bare  wire  can  be  found  from  transmission  line  theory: 

T  (Eq.  27) 

lo  slnh7L 


(1)  Open  circuit 

!(*•)  = 
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(2)  Short  circuit 


T  co*l\L'y(L+£-z')]  (Eq.  28) 

l°  cosK  yL 


The  Herz  potential  of  the  bare  wire  is  found  using  each  of  these  currents.  When 
added  to  the  Herz  potential  of  the  insulated  part  the  following  results  are  obtained: 


U)  T-r^)  r,  l  e~yf!  (.  ccJSXL  - 1  <Eq.  29) 

4-jra-R  l  'U  sink  yL ) 


(2)  ttxW  -  X  l  g-1rR  ( 1  +  t^k'YL  <E9-  30> 

1  ~  4-1 Vo-K  l.1  / 

If  HfLi  ^  1  ,  the  term  in  brackets  in  Eq.  29  is  approximately  equal  to 

l+(L/2  £  ).  Similarly  for  Eq.  30  the  term  in  brackets  is  approximately  equal 
to  1  +  L I  £  .  Comparing  these  results  with  Eq.  25,  it  follows  that,  for  the 

short-circuited  coaxial  antenna,  the  equivalent  dipole  moment  is  equal  to: 

(r  l) 5C  =  T^+CXL)  (Eq.  31) 

where  the  factor  depends  on  whether  the  bare  end  must  be  considered  as 
a  short-circuited  transmission  line  ( =  1),  an  open-circuited  tr  ansmission 
line  (  CX  =  0..5),  or  something  in  between  (0.  5  <CO(  l).  In  designing  an  antenna 
one  would  generally  select  the  case  of  £<-=  l  and  determine  L  as  a  function  of  £ 
in  such  a  way  that  the  bare  part  constitutes  a  short-circuit  for  the  coaxial 
antenna.  If,  for  example  L^>  ^^vWoKCf  and  ^=100  m,  then  it  was  found 
that  L  =  25  m  would  satisfy  the  conditions  for  frequencies  of  60  cps  and  above. 
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3.  ESTIMATE  OF  THE  ACCURACY  OF  FIELD -STRENGTH  CALCULATIONS 


AND  MEASUREMENTS 

The  expressions  for  the  field  from  a  horizontal  dipole  in  the  conducting  medium 
have  been  treated  in  section  2. 1.  As  mentioned,  these  expressions  are  valid 
for  a  conducting  half-space,  and  corrections  for  discontinuities  were  derived 
in  section  2.  2.  In  all  cases  the  expressions  are  the  result  of  approximations 
introduced  to  facilitate  the  mathematics  and  to  obtain  solutions  in  a  closed  form. 
From  the  point  of  view  of  the  experimenter,  these  various  approximations 
obscure  the  picture  and  make  it  difficult  to  determine  whether  deviations  from  a 
theoretical  result  are  due  to  measuring  errors  or  to  a  built-in  inaccuracy  in  the 
expressions. 

An  estimate  of  the  accuracy  of  these  various  expressions  has  therefore  been 
made.  Such  an  investigation  is  rather  involved,  and  the  estimate  was  performed 
only  for  the  horizontal  electric  component  in  the  conducting  medium.  The  results 
are  presented  in  Figs.  13  and  14. 

The  ft  -z  half-plane,  where  ft  and  z  are  expressed  in  terms  of  skin-depth,  is 
basically  divided  by  two  curves  A  and  B.  At  points  on  these  curves,  the  modulus 
of  the  field  due  to  direct  transmission  is  equal  to  the  modulus  of  the  field  due  to 
the  transmission  path  along  the  surface  and  perpendicularly  downwards.  The 
actual  field-strength  measured  on  the  curves  depends  on  the  phase  relationship 
between  the  two  transmission  paths.  This  phase  relationship  can  create  maxima 
and  minima,  and  the  dotted  lines  indicate  where  one  or  the  other  of  the  fields  is 
predominant  and  where  the  expressions  shown  can  be  used  subject  to  the 
approximation  introduced.  The  error-limits  have  been  chosen,  rather  arbitrarily, 
at  t  4%  and  "t  40%,  and  thereby  suggest  which  regions  should  be  avoided  during 
experiments  if  one  wants  to  check  measured  results  againpt  theoretical 
expressions. 
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4.  EXPERIMENTS  CONCERNING  EM-FIELDS  AND  PROPAGATION 


Propagation  experiments  and  the  measurements  of  field -strength  have  been 
mentioned  several  times  in  this  report.  There  were  several  reasons  for 
performing  these  experiments.  Firstly,  the  difficult  relationships  between 
the  various  parts  of  the  EM  field  makes  this  an  interesting  subject  for 
experiments  that  provide  a  good  introduction  to  the  whole  topic. 

The  second  reason,  which  may  be  of  even  greater  importance,  is  related  to  the 
measurement  of  EM  background  noise.  If  the  background  noise  is  to  be 
measured  at  absolute  levels,  some  kind  of  calibration  of  the  measuring  equipment 
is  needed.  In  general,  it  is  quite  easy  to  determine  the  amplification  of  the 
equipment  by  introducing  a  known  signal  across  the  detector  and  measuring  the 
output,  but  this  system,  which  is  widely  used  by  experimenters,  does  not  take 
the  detector  itself  into  account.  One  exception  is  the  magnetic  detector,  which 
can  be  surrounded  by  a  coil  or  a  set  of  coils  {for  example,  Helmholz-coil  systems) 
with  which  a  known  field  is  created.  However,  this  is  not  suitable  for  use  at  sea, 
nor  is  it  applicable  to  extended  electrode  systems.  The  propagation  experiment 
is,  therefore,  a  means  of  calibrating  the  measuring  system  on  the  basis  of  a 
known  dipole  moment  and  the  known  relative  positions  of  transmitting  dipole  and 
receiver. 

It  should,  however,  be  kept  in  mind  that  this  way  of  calibrating  is,  in  fact, 
only  good  for  nearby  signals,  and  the  question  arises  as  to  what  conclusions  can 
be  drawn  with  respect  to  the  source  of  background  noise  arising  from  far-off 
signals  on  the  basis  of  the  calibration  with  nearby  signals.  In  clean  geometry, 
the  answer  is  that  it  should  be  fairly  safe  to  calibrate  the  detector  in  this  way 
when  far  from  coastlines  and  over  deep  water.  The  polarization  of  the  back¬ 
ground-noise  field  and  the  orientation  of  the  receiver  mxist,  however,  be  kept  in 
mind. 
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The  accuracy  of  such  a  calibration  depends  on  a  number  of  factors,  one  being 
the  error  of  approximation  in  the  expressions  for  the  field-strength,  as 
already  mentioned  in  Section  4.  Other  factors  deal  with  the  relative  positions 
of  transmitter  and  receiver,  since  the  angle  dependence  of  the  components 
requires  a  knowledge  of  the  relative  angle  between  transmitter  and  receiver, 
as  shown  in  Section  2.1  A  further  source  of  error  is  the  horizontal  distance  /D  . 
As  seen  from  various  expressions,  P  generally  appears  m  terms  of  jO"  , 
where  n  is  usually  3.  This  means  that  a  relative  error  of  10%  in  the  meas¬ 
urement  of  ft  introduces  a  relative  error  of  30%  in  the  field-strength  measured. 
From  a  practical  point  of  view,  this  source  of  error  is  probably  the  mo3t  serious 
one,  and  has  to  be  taken  into  account  when  comparing  experimental  results. 

An  example  of  a  typical  propagation  experiment  is  shown  in  Figs.  15  and  16. 

The  experiment  was  performed  with  a  dipole  excited  by  short  pulses,  as  shown 
in  Fig.  15.  The  pulses  are  generated  in  a  transistorized  pulse -generator 
triggered  by  a  low-frequency  oscillator.  The  width  of  the  pulses  at  their  base  is 
about  40  ms  and  does  not  vary  greatly  when  the  pulse  repetition  frequency  is 
changed.  The  relative  level  of  harmonics,  therefore,  varies  with  frequency. 

In  the  present  case  the  repetition  frequency  was  10  cps.  The  upper  cut-off 
frequency  of  the  equipment  was  set  at  about  30  cps,  so  that  only  the  10  and  20  cps 
components  are  of  importance. 

The  experiments  were  performed  in  the  following  way.  The  receiver  electrode 
system  was  suspended  from  a  raft  to  a  depth  of  about  10  m,  which  corresponds 
to  1/8  and  1/5  of  a  skin-depth  of  10  and  20  cps  respectively.  No  corrections  were 
made  for  the  attenuation  due  to  the  depth.  The  transmitting  dipole  was  towed 
behind  a  small  boat  on  a  constant  course  and  at  a  constant  speed.  An  attempt 
was  made  to  keep  the  course  parallel  to  the  receiving  electrode  system,  the 
distance  of  closest  approach  being  40  m.  The  received  signal  was  recorded  on 
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tape,  and  the  rms -values  of  the  frequency  components  of  the  field -strength 

were  determined  at  different  points  along  a  track  by  power  spectrum  analysis. 

The  rms -values  of  the  dipole  moment  at  the  various  frequencies  were  found 

from  a  spectrum  analysis  of  the  transmitting  current.  The  transfer-function 

'  E/M  I  was  then  computed  in  rms  values  of  —i— ~ — 

1,1  amp-m 

The  experimental  points  are  shown  in  Fig.  16,  together  with  the  curves 
representing  the  transfer -function  under  the  given  geometrical  conditions. 

There  is  reasonable  agreement  between  experiment  and  calculation.  It  should, 
however,  be  kept  in  mind  that  in  this  case  the  calculations  have  been  performed 
for  a  conducting  half-space.  In  fact,  the  water-depth  at  the  measuring  place 
was  only  100-150  m.  This  fact  could  eventually  explain  why  the  20  cps  points 
fit  the  curve  somewhat  better  than  the  10  cps  points. 

The  propagation  experiments  described  made  use  of  the  latest  equipment  developed 
by  the  EM  Group.  This  equipment  will  be  described  in  detail  in  the  second  report 
of  this  same  series  (T.R.  37). 
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Fig.  1 

Skin-depth,  S  ;as  a  function  of  the  relevant  frequencies  for  various  conductivity  values  typical 

of  sea-water 
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Geometry  of  the  measurements  of  horizontal  field  strength;  seen 


from  above  the  surface 
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Fig.  6  The  angle  dependence  of  the  horizontal  E  components  at 
the  surface  and  where  I'f/o  |»  l- 
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The  influence  of  the  sea-bed  on  propagation 
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Fig.  8 


Correction  factor  for  use  with  limited  water-depths 
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An  infinitely  thin  conducting  plate  with  a  round  hole  —  used  as  a 
model  of  an  insulating  island  in  an  ocean 
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West-east  cross  section  of  the  island  of  Tino  (Gulf  of 
La  Spezia),  showing  positions  of  the  two  stations  and 
the  approximation  used 
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Values  of  H  / H  and  E  /E  at  various  points  on  the  model 
x  o  z  o 

of  the  island  of  Tlno 
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Estimate  of  errors  in  field  expressions  for  E 1  p  of  dipole  at  the 
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Waveform  of  transmitter  current 
Peak  value  ~  10  amp-m 
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Fig.  16  Result  of  propagation  experiment. 


APPENDIX  I 


Computer  programme  for  evaluating  the  modulus  of  the  transfer- 
function  in  the  range  of  -200  —  x^  —  200  m  for  various  combinations 
of  the  frequency  f,  the  distance  of  closest  approach  t  the 

conductivity  and  the  angle  0  .  (Written  in  "Algol"  for  an  Elliott 
503  computer):- 


Transfer  function,  modulus  of  E/M,  in  very  near  field  close  to  boundary 
of  lossy  half  space,  issue  No.  2,  L.B. -N. ,  Nov.  1964; 
begin  real  f,  rhoo,  sigma,  theta,  xl,  t,  u,  v,  P,  Q,  rho,  TF;  integer  c; 
switch  s:  =  next; 

procedure  PQ  (x,  y,  z,  P,  Q);  value  x,  y,  z;  real  x,  y,  z,  P,  Q; 
begin  real  p,  q,  r;r:‘=:0.  0019869*x*sqrt  (y*z); 

p:=:exp(-r)*sqrt(2*rT  2+2*r+l); 
q:=arctan  (r/(l+r))-r; 

P;=p*cos(q);  Q:=p*sin(q) 
end  of  procedure  PQ; 

next:  read  f,  rhoo,  sigma,  theta;  samelme;  print ££15 ?f  =  ?,aligned  (2,  1),  f, 
£cpsrhoo=?,  aligned  (3,  1),  rhoo,  £,m  sigma  =?,  aligned  (1,  1),  sigma, 
5,S/m  theta  =?,  aligned  (2,  1),  theta,  Jc  deg.  ?;  print  £.5-1??; 
t:=theta*0.  0174533;  u:=cos(t);  c:  =  0;  prefix(jL3t  s3?  ? );  scaled  (3); 
for  xl;  =  200  step  10.  0  until  200.  5  do 

begin  if_c-(c  div  5)*5  =0  then  print  1  ?  ?  ; 

rho:  =  sqrt  (xlf"  2+rhoo^2);  v:  =*arctan(xl/rhoo); 
c;=c+l;  PQ(rho,  f,  sigma,  P,  Q); 

TF;=sqrt  ( (2*P*u-u-3*cos(v-t)t  2+(2*Q*u)f  2)1  rhof  3/ 

12.  56637/sigma; 
print  TF 

end; 

goto  next 

end  of  program  transferfunction; 
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